D mesons in matter and the in-medium 
properties of charmonium 



Bengt Friman", Su Houng Lee"'^'^ and Taesoo Song^ 

" GSI, Planckstr. 1, D- 64-291, Darmstadt, Germany 
^ Department of Physics and Institute of Physics and Applied Physics, Yonsei 

University, Seoul 120-749, Korea 
Cyclotron Institute, Texas A&M University, College Station, Texas 11843, USA 



Abstract 

We study the changes in the partial decay widths of excited charmonium states 
into DD, when the D meson mass decreases in nuclear matter, taking the internal 
structure of the hadrons into account. Calculations within the 3P0 model for -0(3686) 
and ^^(3770) imply that naive estimates of the in- medium widths based only on 
phase space are grossly exaggerated. Due to nodes in the wave functions, these 
states may even become narrow at high densities, if the D meson mass is decreased 
by about 200 MeV. For the x states, we generally expect stronger modifications 
of the widths. The relevance of the x widths for J/^lJ suppression in heavy ion 
collision is discussed. These phenomena could be explored in experiments at the 
future accelerator facility at GSI. 
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1 Introduction 



The study of hadrons containing heavy quarks is an important tool for un- 
raveling the properties of the non-perturbative QCD vacuum. In particular, 
in mesons consisting of one heavy and one light quark, such as the D meson, 
the heavy quark acts as an almost static source for the light quark, which in 
turn probes the QCD vacuum. This is the basic constituent quark picture of 
heavy-light meson systems within the formulation based on heavy quark sym- 
metry [1]. Therefore, changes in the vacuum condensates, e.g., at finite tem- 
perature and density, are expected to affect the light quark and consequently 
also to modify the in-medium properties of the D meson. This is explicitly 
borne out in the model calculations based on QCD sum rule analysis [2,3] and 
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the Quark Meson coupling model [4], both of which show a reduction of the D 
meson mass of about 50 MeV at nuclear matter density. Lattice gauge theory 
results for the heavy quark potential at finite temperature, suggests a similar 
drop of the D meson mass at finite temperature [5] . 

The reduction of the D-meson mass obviously has important phenomeno- 
logical consequences, ranging from the possible existence of charmed mesic 
nuclei [4] to enhanced subthreshold production of open charm in pA reac- 
tions [6] . Moreover, a reduction of the D-meson mass in matter has also direct 
consequences for the production of open charm [7] and the suppression of 
J/'?/'- mesons [8] in relativistic heavy ion collisions. Such phenomena could be 
explored in experiments e.g. at the future accelerator facility at GSI [9]. 

The decrease of the D-meson mass in matter is by and large due to the restora- 
tion of chiral symmetry in the nuclear medium. This is demonstrated by QCD 
sum rule calculations of the D meson mass in vacuum [10] and in matter [2,3]. 
In the D-meson sum rules, the leading term is dominated by chiral symme- 
try breaking operators. This is contrary to heavy quarkonia, such as J/ip and 
t/'(3686), where the constituent quark and antiquark probe only the gluonic 
condensates. Consequently, such states are expected to experience only a rel- 
atively small mass shift in medium [11-16]. This implies that, with increasing 
density, the DD threshold [2,5,17,18] will cross the energy of some of the ex- 
cited charmonium states. Consequently, the reduction of the DD threshold in 
matter may be reflected in the dilepton spectrum oi pA or AA reactions as an 
increased width of the peaks corresponding to the ■0(3686) state. Furthermore, 
since more than 40% of the J/ ip produced in heavy ion collision emanate from 
the ■?/'(3686) and x states [19], such crossings will induce a stepwise suppres- 
sion of J/ tjj signal due to the successive melting of excited charmonium states 
at finite density and temperature [2,5,17,18]. 

In this letter, we point out that the level crossing between the charmonium 
and the DD threshold does not result in an immediate dissolving of the char- 
monium states, as found in naive calculations, where the participating mesons 
are effectively treated as point particles. When the internal structure is taken 
into account, the effective coupling to the DD final state depends strongly 
on the Q-value, and consequently on the momentum carried by D meson in 
the charmonium rest frame. The overlap of the D mesons with the wavefunc- 
tion of the initial heavy quarks of the charmonium depends strongly on the 
relative momentum of the D mesons in the final state. In other words, the 
effective ipDD coupling constant is sensitive to the wave functions and the 
momenta of the particles involved in the decay. This mechanism provides a 
viable interpretation [20] of the experimentally observed branching ratios for 
the decay of the ip{AOAO) into DD,DD*,D*D*. There are strong deviations 
from predictions based on naive quark spin counting. It was shown that the 
effective coupling between -0(4040) and the D mesons varies rapidly with the 
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relative momentum in the final state. The matrix elements may even vanish 
at certain momenta, corresponding to nodes in the wave function [20]. This 
mechanism was confirmed also in a more sophisticated potential model for the 
heavy quarkonium [21]. 

As we will show, a similar mechanism is active also in the decay of iIj{368Q), 
'0(3770), Xeo(3417) and Xc2(3556) into DD as the mass of the D meson de- 
creases. Due to the nodes in the radial (■^(3686)) and orbital (■0(3770)) wave 
function, the partial decay widths first increase and then decrease as the mass 
of the D meson is reduced. In particular we find that, for a mass shift of 
200 — 250 MeV, the branching ratio into the DD channel vanishes and then 
increases again when the D mass is reduced further. The resulting widths are 
much smaller than those obtained in the naive picture, where the width is 
enhanced due to the increase in phase space, while the coupling constant is 
kept fixed. 

For the x mesons the picture is somewhat different. The partial width of 
the Xco(3417), increases very rapidly as the DD channel opens up, and then 
approaches zero as the mass of the D meson is decreased even further. On 
the other hand, for the Xc2(3556) the partial decay width increases slowly and 
monotonically, because there is no node in the wave function. 



2 Charmonium states 

We use the harmonic oscillator potential to model the bound state wave func- 
tions and the 3P0 model to describe the charmonium decays. In this ex- 
ploratory calculation our aim is to determine the in-medium properties of 
the excited charmonium states on a semi-quantitative level. Our main result, 
the medium dependence of the effective coupling constants, is due mainly to 
the node structure of the wave functions and does not depend strongly on 
the details of the model. Hence, a more sophisticated calculation based on 
a refined potential will not change the main conclusions of our work. In the 
harmonic oscillator potential model, the wave function of a heavy quarkonium 
state is of the form. 



(f)N,i = Normalization x Y;"^(^, 0)(/3V)^'e-^^''''L5Jii(/3V^) (1) 

where [3"^ = Muo/Ti characterizes the strength of the harmonic potential, M = 
^rric is the reduced mass of the charm quark anti-quark system, and L^^z) is 
a Laguerre Polynomial. The energy eigenvalues are 
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Table 1 

Charmonium states in the harmonic oscillator model. 
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meson A 




meson C meson B 

Fig. 1. Diagrammatic representation of meson A decaying into mesons B and C by 
creating a 3P0 quark anti-quark pair. 

En = huj{n+^), n = 2k + l + l = 2{N + l + (2) 

where, is the number or nodes in the radial direction, including one at infin- 
ity. In table I, we summarize the quark-model assignments of quantum num- 
bers and other relevant information on the lowest lying charmonium states. 
For the ip{3770), the coupling constant g^DD is determined by fitting its width 
using Trot = T^^dd = (4^^/2471) ((m^ - Amlf/yml). For the J/^(3097), 
which is far below the DD threshold, we use vector meson dominance for 
the electromagnetic current of the D-meson, assuming that the form factor is 
dominated by the J/ip at small g^. In vacuum the DD threshold is located at 
3.74 GeV. 



3 The 3P0 model 

We compute the decay of the charmonium states into DD mesons using the 
3P0 model [22]. In this model, the decay of a meson A into mesons B and C 
involves the following invariant matrix element 
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MA-.BC^{A\l[qsqs] '^°\BC), 



(3) 



where 7 is a coupling strength, which characterizes the probabihty for creating 
a quark anti-quark pair in the 3Po state. In fig. 1 we show a diagrammatic 
representation in the rest frame of the meson A consisting of a heavy quark 
anti-quark pair with momenta kq — ks and ks — kq. The decay products are the 
mesons B and C with total momenta /c^ and —ks respectively. The meson B 
is composed of a heavy anti-quark (momentum ks — k^) emanating from the 
parent meson A and a quark (momentum kg) from the 3P0 pair. Now, since 
the decaying quarks has to form the physical meson, the matrix element will 
involve the following overlap integral. 



Ma^bc oc J d%(t)A{'^kq - 2kB)(pB{'2kq - kB)(t>c{'^kq - ke) 

X [uk„sv-kj'''° (4) 

where the mesonic wave functions are the Fourier transforms of the spatial 
wave functions in Eq. (1). The momentum space wave functions are poly- 
nomial functions in the relative momentum multiplied by gaussians. Finally, 
[ukq,sV-kq,s\ OC kqY{"{kq) is the wave function of the produced qq pair. 

The overlap integral in eq.(4), is a function of the meson momentum ks- When 
there are nodes in one of the wave functions (f>, the integral can vanish for 
certain values of ks- This is seen by combining the three gaussians in eq. (4) 
into one. To this end one makes a change of variables 



, 1 + 

kq = kq- ^ _^ 2r2^^' 

where r = |, 5 is the strength of the HO potential of the parent meson A 
and P is that of the emitted mesons B and C. We introduce two independent 
parameters (3 and S to allow for different sizes of the wave function for the 
initial meson and the outgoing mesons. Now, if one or more of the mesons 
involved is in an excited state, the polynomial in the corresponding wave 
function (p is, after the change of variables (5), a polynomial in ks and k'^. 
Consequently, the matrix element (4) is proportional to a polynomial in ks, 
which may vanish for certain values oiks- 

For the case of interest, where the charmonium decays into two pseudoscalar 
mesons DD, the relevant formulae in the 3P0 model can be found in ref. [23] 
for j3 = 5. Here, we re-derive the formula in 3P0 model, allowing for j3 ^ 5. 
The invariant matrix element for the decay A ^ B -\-C is given by. 
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= :;^yT^'PLs{x,r)e-^^ x + /(rf^)]- (6) 

In our case, the flavor factors h[I{di) + I{d2)] — |. Furthermore, 



x^-x ^ml/A-ml, (7) 

which is the scaled momentum carried by the decaying mesons in the rest 
frame of the parent meson A. The decay rate is then given by 



r{A^B + C)^2nP-^^^\M,s\'. (8) 
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We now present the resulting decay rates for the different charmonium states. 
The quark-model assignments of these states are given in Table 1. For r — 1 
{13 — 5), our results reduce to those given of Barnes et al. [23]. 



;i) X(3417) 



r(x(3417) ^ C + B) = 2»3(-^)'x 



^^^~ 3(l + 2r^) '^^ e-^(^. (10) 



(2) x(3556) 



,(1 3p.-.i5o+^5o) _ 1 rV225(i+r2) 



r 



2\2 



r(x(3556) ^ L> + L>) = T^-^rTT nVo ^{r (12) 
(3) V(3686) 



(2 35,^150+^50) _ 2^/^(1 -3r^)(3 + 2r^) 
3(1 + 2r2)7/2 

2r2(l + r2, 

1 + Tr^^I^TT^^^I^TT^^^' > (13) 



(l + 2r2)(3 + 2r2)(l-3r2)' 
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M^(3686)^ 3^ (l + 2r2)7 



r(^(3686) ^D + D)^ ' -D ^2^ i^^^' ; ) ^3 



(4) V(3770) 



(1 3,,,^.5o+^5o) _ _2VTo /^_y/2 / _ l + r2 

" 3 U + 2r2j "^l^ 5(l + 2r2)'^ ^^^^ 



„l/2 p2 rjllcr . s 7 

r(^(3770) - + = ^^7^ (t^) x3 

^V(3770) 32 V 1 + 2r^ / 



x l-— — -xM e ^0^. (16) 

V 5(l + 2r2) ) ^ ' 

Note that in each case, x is defined by eq.(7) with mj^ the mass of the cor- 
responding charmonium state and = mi). The zero in eq.(14) is due to 
the nodes in the radial wave function, whereas those in eq. (10) and eq. (16) 
results from the orbital part. Nevertheless, for r — 1, the two widths in eq.(14) 
and eq.(16) have the same functional form. 

Our model has three parameters, namely, (3, r = 6/(3 and 7: 

• f3 determines the size of the harmonic oscillator potential of the D mesons, 
while 6 is the corresponding parameter for charmonium. Wc fix /? and the 
ratio r by fitting the partial decay width of ■0(4040) to DD, DD* and 
D*D* to the experimental ratios of 1:20:640, as has been done in ref. [20] 
for r = 1. The formulas of ref. [20] can be generalized to r 7^ 1 by the 
following replacement, 

35 / 4 „ 4 A /15 1 + r2 5r2(4 + r2) 



4 • 3H 15 315 ; V 8 1 + 2r2 (1 + 2r2)3 
r2(5 - 9r2 - lOr^) ^ r%l + r^) 4 



^ 2(l + 2r2)4 ^ + 2(l + 2r2)5^ ) ^^^^ 
The fit yields, 

p = 0.3 GeV, and r = 1.04 . (18) 

The resulting oscillator parameter /3 is close to the value obtained in ref. [20] 
{P — 0.31 GeV) and the ratio r is close to unity. The implied difference in 

size of the D meson and charmonium is very small and, in view of the 
simplicity of the model, presumably insignificant. Within the harmonic os- 
cillator model, one can also relate (3 to the mass splitting between the char- 
monium states if one assumes an effective charm quark mass rric- In this 
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Fig. 2. Partial decay width of charmonium decaying into DD, as a function of the 
D meson mass. 

case, for rric — 1.5 GeV, we find P ~ 0.43 GeV. However, since we are inter- 
ested in tfie partial decay widths of charmonium states, we will choose the 
value given in eq.(18). 
• 7 determines the strength of the 3P0 vertex. This parameter may be deter- 
mined by fitting the r[ijj[3770) — > DD), which implies 7 = 0.281. 

These values for the parameters are qualitatively consistent with those ob- 
tained from a best fit to the decays of excited states of light mesons [23], 
p = 0.36 GeV and 7 = 0.4. 

In fig. 2 we show the dependence of the partial widths for the decay into 
the DD channel on the D meson mass in medium. For the J^*^ = 1 
states, we note that when the D-meson mass is reduced, Tdj^^iP (3770)) and 
Tdq^iP (3686)) first grow to about 90 MeV for a mass shift of about 50 MeV. 
Due to the nodal structure of the wave function, a further reduction of the 
£)-meson mass leads to a strong decrease of the partial width, which then 
vanishes for a mass shift of about 200 MeV and 250 MeV, respectively. In 
both cases, the increase of the width beyond this point is slow, so that even 
when the D meson mass is reduced by 400 MeV, the widths remain below 350 
MeV. In fig. 3, we illustrate the important role of the internal structure of the 
mesons by comparing the in-medium width of the ■0(3770) (eq.(16)) with the 
naive estimate obtained by treating the mesons as point particles 



r,^(^(3770)) = -1^^^^ ■ (19) 

Here the coupling constant is kept fixed at its vacuum value QipoD = 15.4. Con- 
sequently, r^^(t/'(3770)) grows rapidly with dropping D-meson mass, due to 
the strong increase in available phase space. Clearly, the width of the ■0(3770) 
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Fig. 3. The solid line is the partial decay width r(^(3770) D + D) from 3P0 
model. The dashed line is assuming a constant coupling normalized to give the 
same leading behaviour at the threshold point. 



in matter is strongly overestimated in the naive calculation, where the internal 
meson structure is ignored. Consequently, the expectation that the states 
would melt instantly as the DD channel opens up is not well founded. 



As shown in fig. 2, the situation for the x states is somewhat different. The 
width of the lightest x meson , the XcO; increases very rapidly as the decay 
into the DD channel becomes possible. Again, a node in the wave function 
leads to a zero in the width as the D-meson mass is decreased further. The 
sudden increase at threshold is due to the fact that in J'^ = 0+, the DD pair 
is in an s-wave state. This implies, that Xco(3417) will dissolve, when the DD 
channel opens up. However, this will probably have only little impact on J/ if) 
suppression in heavy ion collisions, since the feeding of J/ijj^s from the decay 
of the Xco(3417) is expected to be negligible [19]. 



On the other hand, the DD width of the Xc2(3556) is suppressed at first, 
because for = 2+ the D mesons are in an / = 2 state. As the D-meson 
mass is decreased further, this partial width increases monotonically, because 
the wave function of this state has no nodes. Thus, for large mass shifts, its 
width is much larger than for the other charmonium states. This behavior 
could have an important effect on J/ip suppression in heavy ion collision, 
since at sufficiently high densities, the produced Xc2(3556) mesons will melt 
instantly and not contribute to the production of J /ip^s in heavy ion collision. 
This would eliminate more than 10% of the expected J/V' yield [19]. 
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4 Conclusion 



We have shown that the internal structure of the mesons play a crucial role in 
calculations of the in-medium widths of charmonium states. In particular, we 
find that the dropping of the D meson mass in nuclear matter does not lead to 
the extreme increase of the decay width of charmonium states, expected from 
naive phase space arguments. In fact, for both t/'(3686) and t/'(3770) we find 
that the partial decay width vanishes for a particular D meson mass shift of 
about 200 to 250 MeV, due to cancellations in the matrix elements caused by 
the nodal structure of the quark wave functions. For a mass shift of this size, 
the increase in phase space, would give r^^(-?/'(3770)) ~ 1.2 GeV for a constant 
matrix element. The enhanced widths of the '^(3686) and t/'(3770) may be 
observable in the spectrum of dileptons produced in heavy ion reactions. As 
can be seen in fig. 2, the partial width for decay into the DD channel first 
increases to about 90 MeV and then drops to zero for a D meson mass shift 
of about 200 MeV. 

In matter the total width of charmonium states is further enhanced by reac- 
tions with the surrounding particles, like e.g. + N ^ Kc + D. The threshold 
for this process is fairly low; it is energetically allowed for charmonium masses 
> 3.2 GeV. The cross section, computed in the quark exchange model, 
is given in ref. [24]. For the -0(3686) the asymptotic value was found to be 
(Tabs (■0(3686) + A/") ~ 6 mb. The relatively small cross section is related to the 
compact size of the charmonium states. The contribution of this reaction to 
the in-medium width of the -0(3686) is, to lowest order in density, given by 



Tat = - = ((T„„,(^(3686) + N)vreiPn). (20) 

T 

where pn is the density of nucleons in nuclear matter and Vrei is the aver- 
age relative velocity in the initial state. Assuming that the -0 has a rela- 
tive momentum with respect to the rest frame of the medium , v^ei = 
| — + where vp is the Fermi momentum of nuclear matter. When 

= 0, this gives the additional width Fat ~ 4 MeV at Pn = Po = 0.17fm~^ 
and F^f ~ 25 MeV at p„ = 4xpo- Even forp^ = 1 GeV/c, we find an additional 
width of only ~ 6.9 MeV at p„ = po a-nd ~ 34 MeV at p„ = 4 x po- 

Although there are no calculations of the corresponding reaction for the t/'(3770), 
one expects a cross section of the same magnitude because the root-mean- 
square radii of the two charmonium states are similar [21]. In QCD [25,26], 
the leading order result for aabs is proportional to (r^). Therefore, the addi- 
tional width of the charmonium states -0(3686) and -0(3770) due to scattering 
off the nuclear medium is expected to be small. Thus, we conclude that the 
total width of these resonances in nuclear matter is dominated by the decay 
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into DD mesons. The small scattering contribution may play a significant role 
only at those densities, where the DD channel is quenched due to the nodal 
structure of the wave function. 

For the x(3556) a somewhat different picture emerges. We find that its width 
increases monotonically as the D meson mass decreases. In hadronic colli- 
sions [19] a large fraction of the J/ip^s stem from the radiative decay of x's. 
If the DD width of the x(3556) increases in matter, the probability for the 
decay into J/ipj decreases. Thus, part of the suppression of J/ip^s in heavy-ion 
collisions may be due to this effect. The increase in the width of the x(3556) 
could be observed by measuring the J/V'T decay [19] in future heavy ion ex- 
periments [9]. 

In addition to the medium effects due to the partial restoration of chiral sym- 
metry, one also expects mass changes due to modifications of the confining 
potential in nuclear matter. In our model, this effect will be reffectcd in a 
density dependence of the oscillator parameters /5 and 6, which parameterize 
the strength of the harmonic oscillator potential. A rough estimate of the ex- 
pected change in 5 at nuclear matter density, can be obtained by relating the 
expected mass shift of the J/i/j to a change in the confining potential. We use 
the relation Smj/^ = ^j^Sp^, which applies to the harmonic oscillator model. 
All model calculations [11-16] yield 5m j/^ ~ -7 MeV at nuclear matter den- 
sity. This translates into Sp — — 5.6MeV, which can be safely neglected, since 
it corresponds to a relative shift in /3 of less than 2 % (see eq. (18)) and an 
even smaller change in the final result in fig. 2. 

We conclude that measurements of the in-medium properties of charmonium 
states can provide valuable information on the characteristics of QCD in dense 
and hot matter. 
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